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Abstract. We investigate the location of the (unique) hot spot in a convex 
heat conductor v^ith unitary initial temperature and with boundary grounded 
at zero temperature. We present two methods to locate the hot spot: the for- 
mer is based on ideas related to the Alexandrov-Bakelmann-Pucci maximum 
principle and Monge- Ampere equations; the latter relies on Alexandrov's re- 
flection principle. We then show how such a problem can be simplified in case 
the conductor is a polyhedron. Finally, we present some numerical computa- 
tions. 

1. Introduction 

Consider a heat conductor having (positive) constant initial temperature while 
its boundary is constantly kept at zero temperature. This physical situation can 
be described by the following initial-boundary value problem for heat equation: 

Ut — Au in 51 X (0, oo), 

(1.1) w = l on 17 X {0}, 

u — on dfl X (0, oo). 

Here Q — the heat conductor — is a bounded domain in the Euclidean space 
M.^ , N > 2, with Lipschitz boundary and u = u{x,t) denotes the normalized 
temperature of the conductor at a point x G ft and time t > 0. 

A hot spot x{t) is a point at which the temperature u attains its maximum at 
each given time t, that is such that 

u{x{t),t) = Taaxu{y,t). 

yen 

If is convex (in this case fl is said a convex body), it is well-known by a result 
of Brascamp and Lieb [T| that \ogu{x,t) is concave in x for every t > and this, 
together with the analyticity of m in x, implies that for every t > there is a unique 
point x{t) e at which the gradient Vu of u vanishes (see also [T3]). 

The aim of this paper is to give quantitative information on the hot spot's 
location in a convex body. 

A description of the evolution with time of the hot spot can be found in jl6j : 
we summarize it here for the reader's convenience. A classical result of Varadhan's 
[20] tells us where x{t) is located for small times: since 

—At log{l — u{x, t)} -> dist(a;, dfl)^ uniformly for a; G fi as t — >■ 0+ 

(here dist(a:,9r2) is the distance of x from dH.), we have that 

dist{x{t),M) ast^O+, 

where 

(1.2) M^{xen: dist(a;, dn) = rn} 
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and 

ro — ma.x{dist{y , dfl) : y £ ft} 
is the inradius of fl. In particular, we have that 

(1.3) dist{x{t),dn) rn as t ^ 0+, 

For large times instead, we know that x{t) must be close to the maximum point 
a^oo of the first Dirichlet eigenfunction ipi of — A. Indeed, denoting with Ai = Ai(ri) 
the eigenvalue corresponding to (pi, we have that e^^^u{-,t) converges to locally 
in as t goes to oo; therefore (see [16]) 

(1.4) x(t) — > Xoo as t -> oo. 



While it is relatively easy to locate the set Ai by geometrical means, (1.4 1 does 
not give much information: locating either x{t) or Xoo has more or less the same 
difficulty. In this paper, we shall develop geometrical means to estimate the location 
of x{t) (or Xoo), based on two kinds of arguments. 

The former is somehow reminiscent of the proof of the maximum principle of 
Alexandrov, Bakelmann and Pucci and of some ideas contained in [19j, concerning 
properties of solutions of the Monge-Ampere equation. The estimates obtained in 
this way are applicable to any open bounded set, not necessarily convex. 

Let 17 be a bounded open set and denote by JC the closure of its convex hull; we 



shall prove the following inequality (see Theorem 2.7) 



(1.5) dist(xoo,a^)>cjv diam(l]) 



[diam(rj)2 Ai(r2)]^ ' 

Here, diam(r2) is the diameter of and cjv is a constant, depending only on the 
dimension N, for which we will give the precise expression; observe that the quantity 
diam(ri)2 Ai(r2) is scale invariant. 

When is convex, more explicit bounds can be derived; for instance, the follow- 
ing one relates the distance of Xoo from dQ to the inradius and the diameter: 



(1.6) dist(a:;oo, 9ri) > Cn tq 



diam(r2) J 



where again Cn is a constant depending only on N (see Theorem 2.8 for its ex- 
pression). We point out that the so called Santald point of n always satisfies (1.61 
hence this can also be used to locate such a point (see Section [2] and Remark 2.121 
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The latter argument relies instead on the following idea from [31 [S] . Let 'W^ 
be the unit sphere in M^. For uj E §^~^ and A € M define the hyperplane 

(1.7) 7r(A, uj) = {xeR^ -.x-Lo^X}, 
and the two half-spaces 

(1.8) TT+{X,iu) ^ {x eM.^ : X -UJ > X} n' {X,uj) ^ {x e : x ■ u < X} 

(here the symbol • denotes the usual scalar product in M^). Suppose 7r(A,a;) has 
non-empty intersection with the interior of the conductor O and set 

n+^^nmT+ix,uj). 

Then if the reflection Tx.io{fl^ ^) of with respect to the hyperplane 7r(A,a;) 

lies in Q, then 7r(A,w) cannot contain any critical point of u. This is a simple 
consequence of Alexandrov' s reflection principle based on Hopf's boundary point 
lemma (see Section [s] for details). 

Based on this remark, for a convex body K. we can define a (convex) set ^(/C) 
— the heart of /C — such that x{t) e <^(/C) for every i > (in fact, we will prove 
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that K, \ <v^(/C) cannot contain the hot spot x{t) for any t > 0). The heart <v?(/C) of 
K. is easily obtained as the set 

^{IC)^f]{7r-{X,u^):Tx.Ul^tJcIC}. 

As we shall see, the two methods have their advantages and drawbacks, but they 



are, in a sense, complementary. On the one hand, while inequalities (1.5) and (1.6) 
are quite rough in the case in which 17 has some simmetry (e.g. they do not allow 
to precisely locate even when 17 is a ball), by the second argument, the problem 
of locating x{t) is quite trivial; on the other hand, while in some cases (e.g. when 



n has no simmetries or dfl contains some flat parts, as Example 4.3 explains), we 



cannot exclude that the heart of /C extends up to the boundary dIC of /C, estimates 



(1.5) and (1.6 1 turn out to be useful to quantitatively bound x{t) away from dJC. 
Thus, we believe that a joint use of both of them provides a very useful method to 
locate x{t) or x^o- 

Studies on the problem of locating a;oo can also be found in [B] : there, by argu- 
ments different from ours and for the two-dimensional case, the location of Xca is 
estimated within a distance comparable to the inradius, uniformly for arbitrarily 
large diameter. 

In Section |4] we shall relate ^(/C) to a function TZ/c of the direction w — the 
maximal folding function — and we will construct ways to characterize it. We will 
also connect TZjc to the Fourier transform of the characteristic function of JC: this 
should have some interest from a numerical point of view. Finally, in Section [5] we 
will present an algorithm to compute TZic when /C is a polyhedron: based on this 
algorithm, we shall present some numerical computations. 



2. Hot SPOTS and polar sets 

In this section, if not otherwise specified, H is a bounded open set and we denote 
by /C the closure of the convex hull of 17. Notice that /C is a convex body, that is a 
compact convex set, with non empty interior. In what follows, \E\ denotes the A'^- 
dimensional Lebesgue measure of a set C and \dE\ the {N - l)-dimensional 
Hausdorff measure of its boundary; also, uj^ will be the volume of the unit ball in 



2.1. Preliminaries. We recall here some notations from '18]. The gauge function 
jp of JC centered at a point p S /C is the function defined by 



jp{x) — min{A > : x — p <E A(— p + /C)}, x G 



Observe that we have jp(t{x — p) + p) — t jp{x) for every i > 0; in particular, if 
e /C then jo is 1-homogeneous. We set 



(2.1) 5p(a 



]p {x) - 1, if x e \ {p} 
— 1, ii X = p. 



so that gp is the convex function whose graph is the cone projecting dK, from the 
point (p,-l) e 

It is also useful to recall the definition of the support function hjc of /C, that is 



(2.2) /ik;(0 = max{a; : X e /C}, 



As it is easily seen, h/c is a 1-homogeneous convex function; viceversa, to any convex 
1-homogeneous function /i, it corresponds exactly one convex body whose support 
function is h (refer to for instance). 
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The polar set of K, with respect to p is the convex set K.* coinciding with the unit 
baU of the "norm'j^ || • ||* ~ hic{-) centered at p, that is 



/C; {y e : (x-p)- {y~p) <1 for every x e IC}; 

Hp is in the interior of /C, then IC* is compact. Observe that this can be equivalently 
defined as 

JCp^lyeR'^ : {x - p) ■ {y - p) < jp{x) for every 2; e M^}. 

We also recaU that for every convex body K, the function i/j : /C — ?> [0, 00] defined 
by ip{x) = |/C*| attains a positive minimum at some point sk: € IC, which is cahed 
the Santald point of K. (see [IS]). When referring to the polar of JC with respect to 
its Santalo point, we simply write /C*, instead of 1C%^: we will see that the method 
developed in the next subsection for estimating the hot spot, applies to the Santalo 
point as well. 

It is not difficult to see that /C* coincides (up to a translation) with the subdif- 
ferential dgp of gp at the point p, i.e. 

(2.3) dgpip) ^ ic; ~ p, 

where for every xq 

dgpixo) = {e G : gpix) > gp{xo) + ^ ■ (x - xo), x G M~}. 

Finally, we will need the following monotonicity property of the subdifferential 
of a fmiction: the proof can be found in [S]. 

Lemma 2.1. Let Ui and U2 be continuous convex functions on K, such that Ui = U2 
on die. Define: 

du,{IC)^ [j du,{x), 1 = 1,2. 
If Ui < U2 in JC, then du2{IC) C dui(K). 

2.2. The polar set of JC with respect to the hot spot. The following result 
holds for a general domain £7 and is the cornerstone of our estimates. 



Theorem 2.2. (i) Let u be the solution of the initial-boundary value problem (1.1 1 
and, for every fixed time t G (0, +00), let x(t) £ be a hot spot at time t, that is a 
point where the value 

M(t) = miixu{-,t) 
n 

is attained. 
Then 



(2.4) |/C:(,)|<[iVAf(t)]-^M \ut{x,trdx, 

Jc(t) 

where C{t) is the contact set at time t, i.e. the subset offl inhere —u{-,t) coincides 
with its convex envelope. 

(ii) Let Ai(f2) and (f>i be respectively the first Dirichlet eigenvalue and eigenfunc- 
tion of ~A in Q. Let x^o be a maximum point of (pi in and set Mao = Pii^oo)- 

Then 



(2.5) \JC* 



< 


r Ai(r!) 1 


"I 









4>i{x)^ dx 



where C is the contact set of pi, i.e. the subset of where —pi coincides with its 
convex envel 



"'^Properly speaking this is not a norm, since in general we have |1 — 7^ II 3^11 *• 
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Proof, (i) For t e {0,oo), let C/(*) : K. ^ R denote the convex envelope of — t); 
then, is a continuous convex function in /C, such that U^^^ = on dIC. 
The function G defined by G{x) = M{t) gx{t){x), for x e /C, is such that 

G > in /C and G = J/^*) on a/C, 



and hence dG{lC) C dU^^\lC) by Lemma 2.1 By the rescaling properties of the 
subdifferential and (2.3), we know that 

dG{K.) = M{t) {K.l(,^-x{t)), 

thus, 

|/c:(,)|<M(t)-^|ac/W(/c)|. 

On the other hand, by Sard's Lemma and the formula for change of variables 
(see for instance O Section 1.4.2]), we obtain 

\dU^'\lC)\< [ \detD^uix,t)\ dx, 



JC(t) 

with C{t) = {x € : U^^''{x) = —u{x,t)}. Observe that the contact set is not 
empty, thanks to the fact that x{t) e C{t) and moreover we have \C{t)\ > 0. Now, 
by the arithmetic-geometric mean inequality, we have in C{t) that 

\^u(xA)\ 



\<letD'^u{x,t)\^'^ < 



N 



which yields 



C(t) 



\det D^u{x,t)\ dx < N 



-N 



|AM(x,t)r dx 



Therefore, we finally obtain that 

i/c:(,)i< [TV M(t)]-^ 



C(t) 



C{t) 



\Au{x,t)\^ dx 



and we conclude the proof by simply using the equation Au — ut- 
(ii) The proof runs similarly to case (i). 



□ 



Estimates (2.4 1 and (2.51 are generally difficult to handle. The following weaker 
forms of (2.5) may be more useful. 

Corollary 2.3. Under the same assumptions as in Theorem \2.^ we have: 



(2.6) 
and 
(2.7) 



1^: 



< 



N 



\n\. 



\c\> 



N 



N 



\JC*\; 



we recall that K.* denotes the polar set of K. with respect to the Santald point. 

Remark 2.4. As is well-known, |/C*| can be estimated from below by mAr/|/C|, 
where mjv is a positive constant (see [Ej). Thus, (2.7) becomes 

1 N 



|/C| 



N 

Ai(r!)|/C|2/w 



Remark 2.5. In fS^, p. 223] the following problem is posed: 
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[...] Suppose u > is a solution of 

(2.8) - Au = f{u) in fl, u = on dfl 

in a bounded domain il in E^, say u E C^{il). Is there some e > 
depending only on H. (i.e., independent of / and u) such that u has 
no stationary point in an e-neighbourhood of dfl7 

This is true for iV = 2 in case f{u) > for u > 0, but for iV > 2 
the problem is open. [...] 

Here, we point out that by the same arguments used for Theorem |2.2[ we can 
easily prove the following estimate: 

where L is the Lipschitz constant for / and xq is a point where u achieves its 
maximum Mq. When /(O) = 0, we obtain the inequality 

(2.10) l^-l^iiv) 



— an estimate, similar to (2.6), that can be used to bound dist(xo, dIC) from below 
in a way similar to that of Theorem |2.7| below. 



An interesting instance of (2.91 occurs when / = 1 — in this case u is the 
torsional creep of an infinite bar with cross-section Q; we thus obtain: 



< 



This inequality can also be viewed as an estimate for the maximum Mq in the spirit 
of the Alexandrov-Bakelman-Pucci principle. 

Using the definition of the polar set, it is easy to see that |/C*| goes to oo as the 
point X approaches the boundary. The following lemma gives a quantitative version 
of this fact and helps us to provide explicit estimates of the position of Xoo ■ 

Lemma 2.6. Let p be any point belonging to the interior of IC and define R{p) = 
max{|p— : y£dlC}. Then 

^ ' ' - i?(p)^-idist(p,a/c) ■ 

Proof. Set d = dist(p, 9/C) and R = R{p). Obviously JC is contained in the ball 
B{p^ R) centered at p with radius R and in the halfspace 

if+ = {yeM^ : {y^p).{p-p)<d'} 

supporting JC at any point p such that |p — p| — d. Set E — B{p, R) n H~^; then 
IC C E whence 

ic;de;. 

Now notice that E* is the convex envelope of the union of the ball B{0, R^^) and 
the point q = p + d^^{p — p); its volume is explicitly computed: 



with a — d/R G [0,1]. Thus, (2.11) is readily obtained by observing that the 
function of a into the braces is increasing and hence bounded below by 1. □ 

We are now ready to prove the first quantitative estimate on the location of Xoo '■ 
this will result from a combination of the previous lemma and (2.6). 
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Theorem 2.7. Under the same assumptions of Theorem \2.S\ we have that 

jN^i . diam(ri) 



(2.12) 



dist(xoo,9/C) > N'^'^LON-i 



r2|i/JV diani(17) Ai(f7) 
In particular, the following estimate holds true: 

diam(f2) 



N 



(2.13) dist(xo„ 9/C) > 2^ iV^-i 

UJN 



diani(17)2 Ai(17) 



N ■ 



Proof. Applying Lemma [2l3| with p = x^o and Corollary |2.3| gives: 

N 

\n\ > 



Ai(») 
TV 



ujN-i/N 

i?(xoo)^-idist(xoo,9/C) 



Thus, (2.121 easily follows by observing that diam(/C) — diam(f7) > R(xoo)- 



Finally, using the isodiametric inequality 



diam(r2) 



N 



in conjunction with (2.121, we show the validity of (2.13). 



□ 



Estimates (2.12) and (2.131 involve the first eigenvalue Ai(rj), which in general 
is not easy to compute explicitly; when fl is convex, we can estimate Ai(ri) from 
above by means of basic geometric quantities, thus providing an easily computable 
lower bound on dist(a;oo, ^fi). 

This is the content of the following Theorem, which represents the main contri- 
bution of this section. 



Theorem 2.8. If ft is convex, then 

Ar2JV-l 



(2.14) 



dist(a;oo, 9ri) > r^ 



Ai(Si) 



N 



(- 



N-l 



\diam(ri) J 



where rn is the inradius offl, Ai(i?i) denotes the first Dirichlet eigenvalue of —A 
in the unit ball and IPR(f2) = |9f2||ri|^/^~^ is the isoperimetric ratio of^. 
In particular, the following bound from below on dist(a;oo, holds true 



(2.15) 



dist(a;oo, i9f^) > rn 



Ai(Bi)^ UJN Vdiam(r2)) 



(2^iV) 



Proof. The proof of (|2.14|) readily follows by combining (2.12) to the following 

Ai(Bi) \dn\ 



upper bound on Ai(f2) 
(2.16) 



N rnM ' 

proved in [4] Theorem 2]. Using (2.14) and the two inequalities 

N-l 



\^>uJNr^ and \dn\<NujN 



diam(ri) 



we end up with (2.15). 



□ 



Remark 2.9. Observe that using (2.13) and the inequality 
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which follows from the monotonicity and scaling properties of Ai, we can infer 



dist(a;oo, cTi) > rj2 



MB, 



UJn 



Vdiani(17) / 



2N-1 



thus providing a lower bound which is strictly greater than (2.15), as long as the 
ratio rf2/diani(r2) is strictly smaller than 1/2 and > 3. 



Inequality (2.161 is in fact a corollary of a sharper inequality holding for star- 



shaped sets. We can then give a refinement of (2.14) which holds in this larger 



class: to this end, we borrow some notations from [4 . 

A set Q is said to be strictly starshaped with respect to a point xq G 51 if it is 
starshaped with respect to xq and if its support function centered at xq, i.e. 



max {y — Xq) ■ x, 
yen 



is uniformly positive, that is ini^^QQ ^n,a;o (^) > 0- Let fi be a strictly starshaped 
set with locally Lipschitz boundary, as in [i] we define 



W{n) = inf j / 



1 



d(j{x) : Xq ^ fl 



Ian hn,xo{x) 

where da denotes surface measure on dVL. According to this notation, [U Theorem 
3] states: 

Ai(Bi) w(n) 



(2.17) 



N 



Arguing as in Theorem 2.8 and using (2.17) in place of (2.16) gives the following 
estimate. 

Theorem 2.10. Let D, be a strictly starhaped set with locally Lipschitz boundary 
and denote by K. the closure of the convex hull offl. Then 



(2.18) 



dist(a;oo, 9/C) > 



2N- 



1 , , / 



Ai(Bi)^ Vdiam(fl)W^(fl) 



1 



w{n) 



Remark 2.11. We remark that (2.18) is sharper and more general than (2.14) 



and it is at the same time more explicit than (2.12), in the sense that, differently 



from Ai(57), the number W{fl) can be computed directly from the support function 
(which exactly determines a convex set). 

Remark 2.12. It is worth noticing that the Santalo point s/c of /C always satisfies 



(2.6) (as well as (2.4) for every i > 0), then it satisfies all the estimates we proved 



for Xoo in this section. In particular. Theorem 2.8 (or Theorem 2.10) can be used 



as well to estimate the location of the Santalo point of a convex set. 

3. AlEXANDROV'S REFLECTION PRINCIPLE 

In this section, for the reader's convenience, we recall some relevant facts about 
Aleksandrov's symmetry principle, which has been extensively used in many situa- 
tions and with various generalizations (see 0| for a good reference). 



For u! G §^ ^, let 7r(A,CLi), 7r+(A,aj), and tt {X,uj) be the sets defined in (1.7) 
and (1.8). Also, define a linear transformation : — >■ M.^ by the matrix: 



.N 



where 5ij is the Kronecker symbol and the Wi are the components of uj. Then the 
application Ta.o. : ^ defined by 

Tx,u;{x) =A^x + 2\uj, X e M^, 
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represents the reflection with respect to 7r(A,a;). As already mentioned, if is a 
subset of K^, we set ^nn 7r+(A,w). 

Proposition 3.1. Let ft be a bounded domain in with Lipschitz continuous 
boundary dVl and suppose the hyperplane tt{X,uj) defined by (1.7 1 has non-empty 
intersection with Q. Assume that T\^u}{^\ C f2. 

If Q, is not symmetric with respect to ■k{X,uj), then 7r(A,a;) does not contain any 
(spatial) critical point oj the solution u of (1.1). 

Proof. For x G i^^j^ and < > the function 

v{x, t) = u{Tx^uj X, t) — u{x, t) 

is well-defined and is such that 

Wt = Aw in X (0,cxd), 

(3.1) v = Q on r2+^^ X {0}, 

v>Q on dd+^^x {Q,oo). 

Hence w > in f^^^^ x (0, oo), by the strong maximum principle for parabolic 

operators (see [TS]). Since w = on [dVL^^ n 7r(A,cj)) x (0,oo), we obtain that 

^ > on it (w is in fact the interior normal unit vector), by Hopf's boundary 

3v Ou 

lemma for parabolic operators. We conclude by noticing that — — = —2 — — on 
(91]+^n^(A,c^)) X (0,^). □ 

With the same arguments and a little more work, one can extend this result 
to more general situations, involving nonlinearitics both for elliptic and parabolic 
operators. As an example, here we present the following result. 

Proposition 3.2. Let f2 and satisfy the same assumptions as those of Propo- 
sition 



3.1 



in particular suppose that T\,u){^X ^■ 
Let u = u(x) be a solution of class C'^{fl) f] C'^{fl) of the system: 

. > Au + f{u) — and u > in fJ, 

^'^ ' u = on dn, 

where f is a locally Lipschitz continuous function. 

If VL is not symmetric with respect to 7r(A,cj), then ■n{\^uj) does not contain any 
critical point of u. 



Proof. The proof runs similarly to that of Proposition |3.1[ the relevant changes 
follow. The function 

v{x) = u{T\,ujx) - u{x), 
defined for x e ^Ix.u, satisfies the conditions: 

Av + c{x)v — in i^^j^, 
^; > on dn+^, 

where the function c(x), defined by 

r .fjuiTx^^ x))- fjujx)) 
c(x) = I u{rx,ux) - u{x) 

y iov u{T\^uix) = u{x), 

is bounded by the Lipschitz constant of / in the interval [0,maXfYw]. Hence u > 
in i^J^, by the arguments used in [3j. Let c~(2:) — max(— c(a;), 0); then 

Av — c~{x)v<Q and ?;>0 in 



for u{Tx,ujx) ^ u{x), 
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and the strong maximum principle can be applied to obtain that w > in fij^- 
The conclusion then follows as before by Hopf 's boundary lemma. □ 

An immediate consequence of this theorem is the following result. 

Corollary 3.3. Let Q and satisfy the same assumptions as those of Propo- 



sition 3.1 Let u\ he the first (positive) eigenfunction of —A with homogeneous 
Dirichlet boundary conditions. 

If Q, is not symmetric with respect to Tr{X,uj), then 7r(A,a;) does not contain any 
critical point of Ui. 

4. The heart of a convex body 

In what follows, we shall assume that JC C is a convex body, that is a compact 
convex set with non-empty interior. Occasionally, we will suppose that /C c is 
of class C^, i.e. a set whose boundary dJC is an {N — 1) -dimensional submanifold 
of of class C^. 

4.1. The maximal folding function. We are interested in determining the func- 
tion given by 

(4.1) 7^;cH :=niin{A€M : TxA^^lJ'^^}, E 

which will be called the maximal folding function of /C; TZf^ defines in turn a subset 
of /C - the heart of /C - as 

<^(/C) = {a; e /C : x ■ uj < 7^K:(w), for every lo e S^"^}. 

Of course, ^(/C) is a closed convex subset of /C. Observe that TZjc can be bounded 
below and above by means of the support functions of <^(/C) and /C : 

(4.2) ;ip(yc) H <nKiuj)<hK:, i^e . 

The following results motivate our interest on <^(/C) and TZjc- 

Proposition 4.1. Let JC he a convex hody. 

(i) The hot spot x{t) o//C, the point Xao and any limit point of x{t) as t ^ 0^ 
always helong to '^{IC); moreover, x{t) and Xoo must fall in the interior of 
^(/C), whenever this is non-empty. 

(ii) The center of mass of K., 

always helongs to the heart <^(/C) o//C. 

(iii) If K, is strictly convex, the incenter xj^ of K. helongs to ^(/C). 

(iv) Let Xfc — 0. // there exist i (1 < I < N) independent directions wi, . . . , 
such that TZfciujj) ^ 0, j ~ 1, . . . , £, then 

e 

9(/C) C /Cn Pi tt{0,ujj). 

In particular, if i — N, then Z'{1C) reduces to xjc and the hot spot of IC is 
stationary. 

(v) Let 



(4.3) rr = max < min 

then 

<^(/C) CB(aJyc,ryc). 
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Proof. Items (i) and (iv) follow by observing that, for A = TZic(u]), the set K.j^^ 
'T\,Lj{^x i^) is contained in K. and is symmetric with respect to 7r(A,w). Hence, 



X — xk. ■ LO = 



\JC\ 



K\{K+ ^un,^{K+ J) 



[X-yuj] dy 



and the last term is non-negative, vanishing if and only if K. is w-symmetric. 
Items (ii) and (iii) are easy consequences of Proposition 3.1 and Corollary 3.3 
For a fixed 6 G §^~^, let us define 

a{0) = max{i : xk + t6 e ^(^)}, 

which is non negative, thanks to (ii). Then x — xjc + a{9)9 G ^(/C) and 

Xk. ■ + a9 ■ Lu < TZ]c (w) , 

for every uj e §^~^ such that oj ■ 9 > 0. Hence 

. 7^K;(w) - Xk; • w 
a < mm , 



thus taking the maximum as 9 varies on S we obtain (4.3) 



□ 



Informations on convex heat conductors with a stationary hot spot can be found 

in [aiiiiioiiiiiiisiiTi]. 



Remark 4.2. Formula (4.3 1 deserves some comments: observe that for every fixed 
6 G §^^^, the minimum problem inside the braces amounts to finding a direction 
Lu close to 0, so to maximize 9 ■ u), and such that at the same time we can fold /C as 
much as possible, so to minimize the difference TZjci^) ~ xic • w. 

We conclude this subsection by an example that shows how the simultaneous 
application of Proposition |4. 1| and the results of Section 2 substantially benefits the 
problem of locating x^o ■ 



Example 4.3. Let us consider a spherical cap 
B'^ = \ {xi,...,xn) e 



N >, 

i=i J 



with < fi < R. Thanks to the simmetry of it is easily seen that its heart is 
given by 



^iB+)^{{xi,...,XN)e 



Xi ^ ■■ ■ ^ XN-1 = 0, H < XN < {R + /^)/2}, 



which is a vertical segment touching the boundary dB^ at the point (0, . . . , 0, 
In particular, by this method we can not exclude that the hotspot x{t) (or the 
point Xao) is on the boundary. However, we can now use the results of Section [2j 
to further sharpen this estimate on the location of Xao- indeed, applying Theorem 
we get 



dist(xoo,9B+) >iR-fi) 



Ai(Bi) 



N 



ujjv-i ( R- fj. 
R + H 



UJN 



where we used that diam(_B+) = 1\J R? — and = (i? — 
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4.2. Computing IZiq. The following theorem whose proof can be found in [31 The- 
orem 5.7] guarantees that, for a regular set (not necessarily convex), the maximal 
folding function is never trivial. 

Theorem 4.4. Let VI he a hounded open (not necessarily convex) subset of M.^ , 
with boundary dft, and denote by K- the convex hull of fl. 

For every uj G S^~^, there exists £ > such that, for every A in the interval 
[hiciuj) — e, h]c{u])), we have: 

(i) r^A^lJ c n; 

(ii) v{x) ■ UJ > Q, for every x <^ dVLC] 7r+(A,a;). 

Unfortunately, the previous result is just qualitative and does not give any quan- 
titative information about the maximal folding function. Moreover, notice that the 
assumption on cannot be dropped, even in the case of a convex domain: think 



of the spherical cap in Example 4.3 for which we have 7?,g+(— cat) = cat). 



In order to compute TZjc, we need some more definitions. We set 

uj^ = 7r(0,a;) 

and for every y e uj^ we define the segment 

<^uj{y) = {x <E K. : X ^ y + tuj, t G M}. 

Then, we denote by V^j : — > uj-^ the projection operator on uj-^, that is the 
application defined by 

Vu,{x) ^ X - {x ■ uj)u!, xeM.^, 

and, for y in the set 

- the shadow of K. in the direction lj - we define: 

o-Lo{y) = min{i G M : y + G /C} and bu:{y) — max{t G M : y + G K.}. 

We say that a convex body K. is ui-strictly convex if dK. does not contain any 
segment parallel to w. If /C is cj-strictly convex, then for every x = y + tu G dIC 
(equivalently y G dSuj{JC)) such that the normal i/{x) to 9/C at x is orthogonal to 
UJ, the set cri^[x) degenerates to the singleton {x}. 

Remark 4.5. We point out that is a convex function on Suj{JC), while b^ is 
concave; moreover, if we set 

graph+(a^) = {{y,tuj) : y G S^{IC), t > a^{y)}, 
graph" (6,^) = {{y,tu}) : y G S^{K.), t < b^{y)}, 

we have that 

graph+(a„) n graph" = JC 
and, as soon as /C is w-strictly convex, 

graph(a^) U graph(6^) = OK., 
where obviously graph(-) denotes the graph of the relevant functions. 

Theorem 4.6. Let IC CR^ he a convex body. For uj G §^ ^ consider the function 
f : Suj{IC) — > M given by 

(4-4) /^(y) = , y G 5„(/C). 



Th 



en 



(4.5) 7^K;(w) = max /„(y). 
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Proof. Observe that 

^t- = {y + tio:ye SUIC),X < t < bUy)}- 

Let A = TZx,{ijj) \ since ui^j ^J) *- then, for every point y+tui with y e 5^(/C) 
and X < t < bu:{y), we have that Tj^^iy + tuj) e /C; in particular, for t = bu{y), we 
obtain that 2A — b^{y) > a^{y) and hence A > /^(j/). Thus, 

TZtcit^) > max /^^(y). 

If j/o G Si^{IC) maximizes by taking A = fui{yo), we see that T-^^{y + tuj) e /C 
for every y £ Suj{IC) and A < t < bu{y). Therefore, Tji^i^j^) ^ ^^'^ hence 

< fMo)- ' □ 

If we now remember that, for a convex domain /C, the quantity 

is the width of K, in the direction oj, we immediately get a nice consequence of the 
previous Theorem. 

Corollary 4.7. Let K. CR^ be a convex body. Then we have the following estimate 
for the width of^{K,) in the direction uj: 

(4-6) < osc 

yeSuj(iC) 

Proof. We first observe that S-^{JC) = Sui{IC), so that 

f-u,{y) ^ -fu>{y), ye5„(/C), 



and (4.1 ) yields 

(4.7) n^i-oj)^- mm fUy). 

yeSu,{IC) 



Then, from the definition of width, using (4.2), (4.1) and (4.7), we get 

<7^K(w)+7^K(-w) 

yeS^iK.) yeS^(K) y€S^(iC) 

thus concluding the proof. □ 



Example 4.8. In general, inequality (4.6) is strict. For example, in consider 
the ellipse given by 

with < & < a. The function TZx: can be easily computed in this case: for every 
Lu — {lui,uj2) £ we get 

— the set {TZtc{uj)uj : w G is the image of a quadrifolium (a rhodonea with 4 
petals) by the mapping {x,y) H> {x/a,y/b). 

Thus, for example, by choosing the direction lu ~ (l/-\/2, 1/a/2), the right-hand 
side of (4.6 ) equals 

,2 u2 I I 



2x/2 Va2 + 62' 
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while clearly the left-hand side is zero since <^(/C) — {(0, 0)}, due to the simmetries 
of /C. 

This example also highlights the interest of the quantity osc^^^y^) — wiy (^jc), 
which can be seen as a measure of the lack of symmetry of /C in the direction of ui. 



The function f^^ in ( 4.4 ) can be explicitly computed by the use of the Fourier 
transform: this is the content of the next result. 

Theorem 4.9. Let JC be a convex body and for lu G §^^^, let be the function 
defined in (4.4 1. 



(4.8) 



fUV) = r; / N ■ , ' y ^ '^"(^)' 



where X/c denotes the Fourier transform of the characteristic function of K- and d^j 
differentiation in the direction lo. 

Proof. For x ^ K. and ^ S we write x = y + tui and ^ = ry + rw, with y G Suj{IC), 
rj G Lo^ and i, t G M. By Fubini's theorem we compute 



(4.9) 

For r = we then obtain: 
(4.10) X^{ii) 



dx 



XK:iy + tuj) e-'^^dt ) e-'^'^dy 

■OO 



'dt 



Tdy. 



[bUy)~aUy)]e-'y^dy. 



Therefore, by the inversion formula for the Fourier transform, we have: 



(4.11) 



1 



X^{tj) e'y-'^dr, 



(27r)^-i 
By ( |4.9| , we also obtain that 

d 



buj{y) - a^iy), yG5„(/C), 

0, X euj^\S^{JC). 



d^X^iO 



dr 

duXiciri) = ~i 



Xiciv + TUj) 



IS^(K) \Ja^{y) 

buj{yf -auiyf 



te-'*^dt e-'^"dy, 



dy, rjeuj , 



and hence 
(4.12) 



(2^) 



2 

0, 



y G 5„(A:), 
yGw^\5„(/C). 



Formula (4.8) follows from (4.12) and (4.11) at once. 



□ 



Remark 4.10. If /C is a polygon, Xjc can be explicitly computed in terms of the 
vertices of /C. Let /C C M'^ be a (convex) polygon with vertices pi , . . . , p„ ; we assume 
that pi, . . . ,pn are ordered counterclockwise and we set Pn+i = Pi- 

Rewriting X/c as a boundary integral (see [5]) by means of the divergence theo- 
rem, we have that 

'^/c(C) = -T7i2 2^ b+1 - PjK'^j • 
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where 

1 \ Pj+i-Pj 



1 / \Pj+i-Pj\ 



.7 = i, . . . ,n, 



is the exterior normal to the j-th side of K.. Also, d^Xjcir]) is easily computed from 
the previous expression: 



wg'"'" -a)-. I 



4.3. Necessary optimality conditions. We conclude this section by presenting 
some necessary conditions for the optimality of /^j. To this aim, we first state and 
prove an easy technical result for the subdifferential of a function. 

Lemma 4.11. Let O C M*^ he a convex open set. Let (p and tjj be a convex and, 
respectively, a concave function from fl to R. If (p + ip attains its maximum at a 
point yo € fl, then 

(4.13) difiyo) C a(-V)(yo). 

Proof. It is clear that both dip{yo) and 9(— V')(yo) are non-empty. Since yo is a 
maximum point, we get 

<p(yo) + V'(yo) > ^{y) + i>{y) for every y e O, 

and hence 

^{y) - f{yo) + £.■ (y-yo) < -~^{y) + ipiyo) + ^ • (y - yo) 

for every ^ e M^, and y e If ^ € 9(p(yo), then we have ^ e 9(— V')(yo)- D 

As a consequence of the definitions of a^^ and bi^, we have that dai^{yo) U 
d{—b)u{yo) = implies that yo belongs to the boundary of iSaj(/C). 
We are now in a position to state a necessary optimality condition. 

Theorem 4.12. Let K C he a convex body and lo e §^~^. Suppose that f^ 
attains its maximum at a point yo € <Sa,(/C), that is 

^k;(w) = /w(yo)- 

Set A = TZjcic^) and for every xq £ 'PZ^ivo) H dK, denote by Xq its reflection with 
respect to the optimal hyperplane, that is Xq = T\,uiXo. 
Then: 

(i) if Xq ^ Xq, we have 

(4.14) A^{N^{x^)) C N^{xo), 
where 

Nk{x) = e \ {0} : a; • ^ = h^i^} U {0}. 

denotes the normal cone of K, at a point x € dK; 

(ii) if Xq = Xq, there holds 

(4.15) AUN,^{xq)}) C N+{xq), 
where 

^Ki^) = U e N^ix) : ^ w < 0} and N+{x) = € N^{x) : ^ w > 0}. 
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Proof. We can suppose for simplicity that w = cat = (0, . . . , 0, 1). We first suppose 
that ?/o is an interior point. Since by its very definition f^^ is the sum of a convex 
function and a concave one, by Lemma |4.11| 

(4.16) daM ^ d[~b^){yo) . 

Let us now set = yo + ^u{ya) ^] the reflection of in the optimal hyperplane is 
= yo + a^(yo)^^ and 

Nk{xo) = {{V, 1) : ^ e 9(-6„)(yo)} N^ix^) = {(77, -1) : 77 e daUvo)} , 



for graph+(a^) n graph (6^) = JC. Then, since A^uiri, 1) = (77, —1), (4.16) implies 
( pll )- 

If j/o is on the boundary of S^{JC), then 9(— 6tj)(yo) or da^jiyo) may be empty: 
this is clearly the case if graph(a(j) or graph(6t^) have some vertical parts. Observe 
that actually we have the following possibilities: 

(1) di-b^){yo) ^ daM = 0: 

(2) di-b^)iyo)^0. 

If (1) holds, then at every xq G 'P^^ivo) H 9/C, the convex body K. has only 
supporting hyperplanes parallel to e^r: these are invariant with respect to the 
action of Auj, so that their reflections are supporting hyperplanes for JC at Xq and 
formula (4.14) or (4.15) easily follows. 

If (2) holds, we have xq = Xq and let us caU ft = JCl^^ U 7a,w(A^a^(^)- Then 

NKixo) C Nnixo) 
and we have N^{xo) = Nj^{xo), so that 

N^ixo) = Nnixo) \ N+{xo) D N^{xo) \ N+{xo) = A^^(xo). 



By observing that N^^{xo) = A^iN+{xQ)) = A^{N+{xo)), ( |4.15| follows. 



□ 



Corollary 4.13. Under the same notations of Theorem \4-lS\ if dJC admits a 
(unique) unit normal v at the point xq, then it admits a unit normal at the point 
Xq too and 

(4.17) Au^i^ix^) = 1^(2:0) . 

In particular, if Xq — Xq we have i^ixo) G w^. 
Proof. It is sufficient to observe that in this case 

Nk{xo) = {^ ■■ i = tiyixo),t>0}U{0}, 



and hence (4.17) is a consequence of (4.14) or (4.15). 



□ 



Using Theorem |4.12[ we obtain an interesting upper bound on the maximal 
folding function for a strictly convex domain, in terms of its support function. 

Proposition 4.14. If K, is strictly convex, then for every oj G §1^^^ 

'yh^iO + vh^iA^^y 



(4.18) 



TZic{uj) < max 



: e G \ {0} 



Proof. First observe that thanks to the 1— homogeneity of the support function, 
the maximization problem in (4.18) can be equivalently settled in 

The strict convexity of /C impfies that hic G C^{R^ \ {0}) (see [15]). Moreover, 
'^hic{0) = X f or eve ry 9 G Ntc{x), with x G dK,. Thus, with the same notations 
as in Theorem |4.12[ x^ = Vhfc{9) for every 9 G Nk:{x^) and using the condition 
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Figure 1. Here w = (1, 0, . . . , 0), a; = V/ik;(0, V = V/ik;(-4^(0); 
the intersection of the two straight hnes corresponds to TZic{uj), the 
dark dot corresponds to ^{x + y). 



(4.14), we also get Xq — Vhic{A^{6)); on the other hand, a^(yo) = Xq ■ uj and 
^uiyo) — ■ uj, which imphes the following 

Liya) = ^ • ^ • 

Hence we can conclude by simply applying Theorem |4.6| □ 



Remark 4.15. If /C is symmetric with respect to a hyperplane orthogonal to ui, 



then equality holds in (4.18 ) and both quantities equal x/c • Otherwise, in general 
inequality (4.18) is strict as Figure 1 informs us. 

Notice that, following an argument similar to that of the proof of Proposition 
|4.18[ we can in fact give a precise characterization of the maximal folding function 
TZic hi terms of the support function h/c. Precisely the following holds 

(4.19) TIkKi^) = max , 

5es;(w) I 

where 

G \ {0} : VhrziC) = ^hK.{A^C) + ^1UJ for some e M} . 

If K, is not strictly convex (and then Hk: is not C^) the above formula still remains 
valid, up to suitably interpreting the gradient of hx, as the subdifferential dh]c- 

5. Numerical examples 

5.1. The case of convex polyhedrons. If /C is a convex polyhedron, then the 
conclusions of Theorem |4.6| can be improved: roughly speaking, we can discretize 



the optimization problem (4.1), by only visiting the projections of the vertices of 



K. on Lo^ . We begin with the following general result. 

Lemma 5.1. Let A and B he convex sets such that A Q B and let x G dA D dB. 

If dA contains a segment t and x belongs to the relative interior of i, then £ is 
also contained in dB. 

Proof. In other words, if A "touches" B from the interior at x and x is contained in 
the interior of some segment on the boundary of A, then the boundary of B must 
contain all the segment at the same. 
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Figure 2. The two cases (i) and (ii). 



Indeed, let tt be a support hyperplane to _B at x and denote by 7r+ the half-space 
delimited by tt and containing B; then tt is also a support hyperplane to ^ at a; 
and A C TT^. Thus, £ C tt'^ while x E irnr 0; this implies £ C tt, since x is not 
an endpoint of £, and hence £ C dB. □ 

Corollary 5.2. Under the same assumptions and notations of Theorem \4-.lS\ if Xq 
belongs to the relative interior of a segment £ contained in dIC, then T\.uj{£) C dIC 
and 

(5.1) ^(^J)=7^;cH for every y e rU^) . 

Proof. The proof follows from the previous lemma by setting B — K, and A = 
^ n T\.uj{lC) and using the definition of □ 

Theorem 5.3. Let xi,...,Xs G he the vertices of an N -dimensional convex 
polyhedron K. C M^, so that 

/C=|a;eM^ : a; = ^ X,x^, with ^ A, = 1. A^ G [0, 1] I , 



For a fixed uj G S"'^ ^ , let f^ he the function defined by (4.4). 

Then TZjc is the solution of the following discrete optimization problem 

(5.2) 7^/c(w) max{/^(yj) : j = l,...,s}, 

where yj = V^^ixj) is the projection of Xj on S^{IC), for every j = 1, . . . , s. 



Proof. By definition (4.1 ) and Theorem 4.6 we know that the value A = TZ]q{lo) can 

+ 

A, I 



possibly be achieved when the boundary of the reflected cap Tx,ui[K,Xu>) tangent 



to that of K. either 

(i) at a point p ^ 7r(A, u), or 

(ii) at a point q € 7r(A,a;) 

(see Figure 2). Thus, the maximum of /(^ is attained at the projection of either p 
or q on Su:(JC). 

Now, let /C be a convex polyhedron. If p is not a vertex of /C, then p belongs to 
the relative interior of some m-dimensional facet of 9/C, with 1 < m < — 1, and 
hence it belongs to the relative interior of a segment £ with (at least) one end at 
some vertex v of dlC. 

By Corollary |5.2[ we then have: 
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Figure 3. The sets K. and 'v'(A^) in two examples. In the first 
picture observe that, by means of (1.6), we also know that Xoo is 
at a positive (and computable) distance from the boundary of /C. 



In case (ii), the m-dimensional facet F of dK. containing q must be orthogonal 
to the hyperplane 7r(A,a;); however, the same argument used for case (i) can easily 
be worked out in F. □ 

The pictures in Figure 3 show two convex polygons with their relative hearts. 
The hearts have been drawn, by using Matlab, by an algorithm based on Theorem 
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